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Explicit computations of Siegel modular forms of degree two. 



Martin Raum, Nathan C. Ryan, Nils-Peter Skoruppa, and Gonzalo Tornaria 

Abstract 

Unlike classical modular forms, there is currently no general way to implement the computation 
of Siegel modular forms of arbitrary weight, level and character, even in degree two. There is 
however, a way to do it in a unified way. After providing a survey of known computations we 
describe the implementation of a class modeling Siegel modular forms of degree two in Sage. In 
particular, we describe algorithms to compute a variety of rings of Siegel modular forms, many 
of which are implemented in our class. A wide variety of Siegel modular forms (e.g., both vector- 
and scalar-valued) can be modeled via this class and we unify these via a construct we call a 
formal Siegel modular form. We define this notion and discuss it in detail. 



1. Introduction 

Computing niodular forms is of broad interest. For example, the explicit computation of 
elliptic modular forms has led to a number of conjectures that in many cases have become 
theorems. Similar, but neither as comprehensive nor as systematic, computations have been 
carried out with Siegel modular forms as well. A consideration of Siegel niodular forms and 
their Fourier expansions aiming at an implementation should suit two purposes. Firstly, one has 
to obtain the generators of a space or ring of modular forms as easily and quickly as possible. 
Secondly, one wants to construct more complicated modular forms and recognize them as a 
linear or algebraic combination of the generators. 

Computing elliptic modular forms of level 1 is rather straightforward. It is well-known that 
the ring of such forms is a polynomial ring generated by the Eisenstein series E4 and Eq. 
Alternatively, the weight k space is spanned by Ek and all the products EiEk-i for i running 
from 4 to fc — 2 as can be proved by means of methods found in |31j . Computations in level 
bigger than one are typically done via the theory of modular symbols. For instance, this is how 
both MAGMA [6] and Sage |52j compute spaces of elliptic modular forms. 

An elliptic modular form can be thought of as a Siegel modular form of degree 1. The 
theory of the computation of Siegel modular forms of even degree 2 is not as general as that 
for degree 1. There is currently no theory of modular symbols. What there is, though, is an 
explicit description of a number of particular rings of Siegel modular forms. More generally, 
by a theorem of Igusa, we know that every ring associated to any finite index subgroup of the 
full Siegel modular group is generated by polynomials in theta constants (cf. |29[ 130) ). This 
theorem has been exploited by Igusa |27[ I28j to explicitly list the 5 generators of the ring 
of Siegel modular forms of degree 2. A more recent application of this theorem in a number 
of papers [2j 1211 123) by Ibukiyama and his collaborators have characterized rings of Siegel 
modular forms for levels up to 4 as polynomial rings in explicitly described generators. These 
generators are typically described in terms of theta constants. 

Another method to compute spaces of Siegel modular forms (and thereby, possibly rings 
of modular forms) is Poor and Yuen's restriction technique (see for example, f43J). In their 
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papers on this topic they are typicahy trying to compute forms in level 1 and degree larger 
than 2, but in principle their technique should work in degree 2. The potential obstacle when 
computing for level bigger than one is that there are several cusps and one needs to know 
the expansion at each of them. The method proceeds roughly as follows: they define a linear 
map from spaces of degree n Siegel modular forms down to degree 1. They search for linear 
relations among the Fourier coefficients of degree 1 forms and pull them back to linear relations 
among the Fourier coefficients of degree n forms. This allows them to compute the dimensions 
of particular spaces of degree n Siegel modular forms and to write down bases for those spaces. 
To determine these bases, they compute theta series, Ikeda lifts, Eisenstein series, and theta 
constants (see Section 

The above two families of results are both for scalar-valued Siegel modular forms. In some 
cases (e.g., |48j). one can explicitly compute spaces of vector- valued Siegel modular forms. An 
example of such a space is given in terms of the 4 Igusa generators of even weight and an 
explicitly described differential operator (see Section [3.8. ip . We limit our attention to spaces 
of vector-valued Siegel modular forms because the ring of such forms is not finitely generated. 

In addition to providing a useful survey of what is known about the computation of degree 
2 Siegel modular forms (e.g., what rings can be computed, what vector spaces, etc.), we 
describe an implementation [46] of the types of Siegel modular forms often used to describe 
generators. We provide algorithms that are used in our implementation in Sage |52| and that 
are improvements (sometimes minor) over the naive algorithms that follow immediately from 
theorems in the literature. Another feature of our implementation is that it allows for algebraic 
combinations of vector- and scalar-valued Siegel modular forms. In fact, our class was designed 
in such a general way that one could use the same basic class to implement a wide variety of 
automorphic forms. 

Often one is interested in computing Hecke eigenvalues associated to modular forms as 
opposed to computing Fourier coefficients. In our work we are mainly interested in computing 
Fourier coefficients but now mention recent work that allows for the computation of Hecke 
eigenvalues. It is possible to compute Hecke eigenvalues by a purely cohomological approach. 
Until now this has been particularly useful in the case of vector- valued Siegel modular forms 
[151 116] . The trace of the Hecke action on the full cohomology can be determined by counting 
curves modulo p. These calculations are demanding in terms of theory and raw computational 
power. 

The paper is organized as follows. In Section [5] we give the definition of Siegel modular 
forms and provide the definition of formal Siegel modular forms: this is what our Sage package 
actually implements. In Section[31 we do two things when appropriate. First, for a wide variety 
of types of Siegel modular forms of degree 2, we give a survey of known theorems related to 
their computation. Second, we give a description of the algorithm we use to compute Siegel 
modular forms of that type and any theorems related to the algorithm. Then, in Section |4] we 
list the spaces and rings that are described in the literature and that we have computed. The 
final two sections deal with the implementation, its details and the choices we have made. We 
also provide links to data for modular forms in those rings for which we have computed data. 
In Section [S] we describe in broad strokes the class we designed, paying particular attention to 
how we multiply Siegel modular forms of different types and to how we model the precision of 
a given form. In Section [S] we briefly point out how other spaces of automorphic forms can be 
modeled via a class similar to the one we describe. 



2. Siegel Modular Forms 

In this section, we define Siegel modular forms of degree 2. If one were to implement this 
definition directly, one would have a rather inefficient class as it would be essentially a dense 
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power series ring in three variables. What we end up implementing in Sage is a formal Siegel 
modular, a notion we define in this section. Also in this section we show that formal Sicgel 
modular forms do model the Fourier expansions of Siegel modular forms rather well. A formal 
Siegel modular form is in some sense more general than a Siegel modular form: the set of 
Fourier expansions of Siegel modular forms is a subset of the Fourier expansions modeled by 
formal Siegel modular forms. We sometimes take advantage of this fact in the implementations 
of our algorithms. At the same time, however, it respects as much structure as needed to store 
the data associated to the modular form in a compact way. 

2.1. Basic DcEnitions 

We recall the definition of a Siegel modular form of degree two. We will consider the 
symplectic group Sp(2,Q) of degree 2 over Q; i.e., the group of matrices M € M(4,Q) such 
that *-MJM = J. Here J = { ) represents the standard symplectic form. The full Siegel 
modular group Sp(2,Z) = Sp(2, Q) n M(4, Z) will be the central object of study. For any 
r C Sp(2, Q) commensurable to Sp(2, Z), let : T — >• §^ be a character and k, j be nonnegative 
integers. Let 



be the Siegel upper half space of degree 2. Assign to C[X,Y]j, the space of homogeneous 
polynomials of degree j, the GL(2, C)-action 



Now we can define Siegel modular forms; 

Definition 1. A Siegel modular form of degree 2, weight {k,j), and character on F is a 
complex analytic function F : IEII2 ^ 'C[X, Y]j such that 



F{gZ) F{{AZ + B){CZ + D)-^) ^ det(CZ + Df {CZ + D) ■ F{Z) 
for aUg=(^£)eF. 



The space of all such functions is denoted M^j (F, -0), where we suppress if it is the trivial 
character and j if it is 0. If j is positive F is called vector-valued, otherwise it is called scalar- 
valued. We use M,(F) :— 0j.Mfc(F) for the ring of (scalar- valued) Siegel modular forms of 
degree 2 on F. 



Common choices of F and V include F = Fo(^) := ( f^^^^ ^2x2 ) n Sp(2, Z), and V ( c B ) — 



such a character ij: is commonly referred to as a Siegel modular form of level i and character 



2.2. Fourier CoefEcients and Formal Siegel Modular Forms 

Having defined analytic Siegel modular forms of degree 2, we now define a formal Siegel 
modular form. Every analytic Siegel modular form F can be identified by its Fourier expansion. 
Because of a certain invariance (see Proposition l2 . ll for a precise statement) that the coefhcients 
of F possess, the form F can be identified by its Fourier coefficients indexed on GL(2,Z)- 
reduced quadratic forms. So, we define a formal Siegel modular form Cp related to as a 
function from integral binary quadratic forms to the representation space of F. 

We denote by Q the set of all integral binary quadratic forms that are positive semidefinite: 



m2 := {Z = X + iY e M(2, C) : *Z = Z,Y>0} 



{A,p)^A-p:=p{{X,Y)A). 



(2.1) 




Q:={f=[a,b,c\ : 



a,b,ce Z, fe^ - 4ac < 0, a > O}. 
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For / = [a, b, c] we use Mf = ( ''f j . The group GL(2, Z) acts on Q by Ma-s = AMf *A. 

Let be a Siegel modular form on F. It is not a serious restriction to assume, as we shall 
do for the following, that P := ^^^^ ^ n Sp(2,Z) C F, since this can always be achieved by 
passing from F to a Sp(2, Q)-conjugate subgroup. Moreover, we also assume that ij: is trivial 
on P. The invariance of F under the action of P implies that F is periodic and therefore has 
a Fourier expansion of the form 

F{Z)= CF{f)e{aT + bz + CT') 

f=[a,b,c]eQ 

= J2 Cp{f)e{tT{ZMf)) 

f=[a,b,c]eQ 

Here Z := ( ^ ^/ ) (r, t' e Hi and z e C), e{x) = e^'^", and the trace of a matrix A is denoted 
by tr A. 

Let A* \A-^). Assume 

(^Q* ) e F for some A G GL(2,Z). (2.2) 

Then, 

l)Z)^ FiA*ZA-') ^ ((f )) det(A)'= A ■ F{Z). 

In other words, equating the Fourier expansion of the left-hand side to that of the right-hand 
side, we get 

Y^CF{f)e{ir{A*ZA-H-If))^^{{^^ \)) Aet{Af ^ A • C7^^(/) e(tr(ZM/)). (2.3) 
/ / 
Using the fact that the trace is invariant under cyclic permutations of the matrices being 
multiplied, we deduce 

lT{A*ZA-^Mf) = tT{ZA~HlfA*) = tT{ZMA-i.f) 

where the last equality follows from the definition of the GL(2, Z) action on Q. Replacing / by 
A - f in the left-hand side of (|2.3p yields the following: 

Proposition 2.1. Let the notation be as above. Then 

CF{A-f) = ^{{ j;o^)) detiA)^ A ■ CfU)- (2.4) 
for all A in GL(2, Z), such that ( 4,* ^ ) € F. 

This allows us to define a formal Siegel modular form. 

Definition 2. Let i? be a module (or ring) with a left GL(2, Z) action, and let x be a linear 
character of GL(2,Z). A formal Siegel modular form with values in R and (formal) character 
X is a map C : Q ^ R for which 

CiA-f)^x{A)A-C{f) (2.5) 

for ah A e GL(2,Z). 

By Proposition 12. 11 we see that this definition models a Siegel modular form well, with the 
formal character given by 

x(A) = V^((f o))det(A)'=. 
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Moreover, it is enough to store the values of the coefficient map C at the GL(2, Z)-reduced 
elements of Q, since the other coefficients can be computed using (|2.5p . 

We note that there are only 4 possibilities for the character of a formal Siegel modular form: 

X € {1, det, tr, deter} 

where a{A) is the sign of A mod 2 acting as a permutation on the set {(0, 1), (1,0), (1, 1)}. 
In fact, the group of linear characters of GL(2,Z) is non-cyclic of order 4 generated by det 
and a. This follows from the fact that the abelianization of GL(2,Z) is isomorphic to the 
Klein four group. Indeed, GL(2,Z) is generated by the matrices E — , T = I) and 

S= -i).From(5r)3 = a.nd ETE = T-'^ we deduce Tr{S) ^ tt{T)-^ and 7r(T) = 7r(T)-i, 
where n is the canonical map onto the abelianization of GL(2, Z). Therefore, the abelianization 
is generated by tt{E) and tt{T), both of which have order 1 or 2. On the other hand, since we 
have four characters, the order of the abelianization is at least four. This implies our claim. 

For a Siegel modular form of level £, the formal character can only be either 1 or det, 
depending only on the parity of the weight k and the parity of the character ip. 

2.3. Formal Constructions 

In order to highlight the usefulness of thinking of Siegel modular forms as purely formal 
objects, we now describe two analytic constructions from this formal point of view. 



2.3.1. Sicgcl $ operator We now describe the Siegel $ operator formally. The map $ 
takes a degree n Siegel modular form to a degree n — 1 Siegel modular form. In particular, 
cusp forms are in the kernel of $ and Eisenstein series are mapped to Eisenstein series. The 
Siegel $ operator for degree 2 Siegel modular forms has a very simple formula in terms of its 
singular Fourier coefficients: 

<^{F):=J2CF{[0,0,n])q". (2.6) 

n>0 



2.3.2. Hecke Operators In this section we give formulas for the image of a Siegel modular 
form of level N and weight {k,j). We give two formulas. The first is for the operator T{p^) 
when p and N are coprime. The second is for the operator U{p) when p divides the level. 

Let F be a Siegel modular form as above. Let C be its corresponding formal Siegel modular 
form and let C be the formal Siegel modular form that corresponds to the image of F under 
r(/). Then, as found in [ig, 

where R{p^) is a complete set of representatives for SL(2, Z)/rg^''(p''); where, for / — [a,6,c], 
[au,hu, cu] = fu ■— f {{X, YyU); where (io,/3 = ( ^ ) ' ^^'^ where the dot denotes the action 
of do^pU as defined in (|2.1[) . We note that if F is scalar- valued this action is trivial. 

For primes p dividing the level of a scalar- valued Siegel modular form F, there is an operator 
U{p). In [4], Bocherer gives an easy explicit formula for the action of U{p) on a Siegel modular 
form. In particular, let C and C be formal Siegel modular forms corresponding to F and the 
image of F under U{p). Then 

C'[{a,b,c)] = C[{pa,pb,pc)]. 
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3. Recipes and Algorithms 



The generators of rings of Siegel modular forms of degree 2 are often given in one of the 
following ways: a linear combination of theta series, a linear combination of products of theta 
constants, or, if the form is of level 1, it may be given in terms of the Igusa generators. 
Additionally, particular Siegel modular forms of degree 2 can also be written in terms of Maass 
lifts, Eisenstein series, or Borcherds products. In this section, we present the definitions of all 
these terms and algorithms or recipes that can be used to compute them (or, at the very least, 
provide references to such methods). 

3.1. Eisenstein Series 

We define an Eisenstein series of degree 2 by 



where the sum is over nonassociated symmetric coprime matrices C,D £ M2 (Z) (here, two 
integral matrices A and B are called coprime if the matrix products GA and GB are integral 
when and only when G is integral) . Earlier work by Maass i32) and Mizumoto |35l 1361 give 
explicit formulas for the Fourier coefficients of Eisenstein series of level 1. We describe another 
approach here. 

In |14) one can find a formula for the Fourier coefficients of a Siegel Eisenstein series of 
degree 2. The formula is in terms of Cohen's function which we define now. Let 



[L{s,{^))Y.d\f^J^{d){^)d-'cJ^-2s{f|d) {iD = Q,l (mod 4) and ^ 

where D has been written as D = DqP with / G N, Dq is the discriminant of (Q)(\/D) and 
ak{n) ^Y.d\nd''- We define Cohen's function by H{k -1,N)= L^nC^ ~ k). 

Now suppose Cfc is the formal Siegel modular form corresponding to the Eisenstein series 
Ek- Then, according to [141 Cor. 2, p. 80], we have 



We observe that the coefficient of a singular form Ci;([0,0, c]) is C(3 — 2k)ak{c), which is the 
cth coefficient of the Eisenstein series of degree 1 and weight 2fc — 2. We note that one can 
compute Cohen's fimction within Sage. 

We remark that the weight k Eisenstein series can also be computed as the Saito-Kurokawa 
lift of the degree 1, weight 2fc — 2 Eisenstein series. In fact, that's where the formulas in this 
section come from. See Section [3.31 for more details on how to compute Saito-Kurokawa lifts. 

One can also consider Eisenstein series with level. Mizuno [371 Theorem 1] gives a very 
complicated formula for Eisenstein series of degree 2 and odd squarefree level. The formulas 
consists of special values of L-functions, Gauss sums, and the like. We refer the interested 
reader to that paper for details. 

3.2. Klingen-Eisenstein Series 

Given a modular form / of degree ni < n and weight k > n + 7ii + 1 the Klingen-Eisenstein 
series of / is a modular form En.m (/) of degree n and weight k such that the Siegel $ operator 



EkiZ) = det{CZ + D)-'' 

{CD} 




if £> ^ 0, 1 (mod 4) 

if i:> = o 
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maps En^m (/) to /. More precisely, it is given by the summation 

En,nAf)^ E /l^' (3.1) 

J>-ferS,"'_^\r(") 

where / is considered a function on ]HI„ -s-^ H„j and T^^^oo is the ni-parabolic subgroup 

ri"U {( c B ) e r(") : = O for * > m, all J, Aj - for i > m, i < m}, (3.2) 

where the zjth entry of a matrix M is denoted My . 

We restrict to n = 2. Then i?2.o(/) = -£-2 and i?2,2(/) = / as can be seen from the easy way in 
which the corresponding ni-parabolic subgroups can be described. So we have to give formulas 
for E2.1 and this was done by Mizumoto in |35l I36j . We can restrict to the case that / is 
an eigenform. Then, we can assume that / is cuspidal, since otherwise i?2,i(/) is a classical 
Eisenstein series. 

Mizumoto gives an explicit formula for the primitive coefficients of a Klingen-Eisenstein 
series in terms of the coefficients of / and special values of three different L-functions. We do 
not present this explicit formula, but note, instead, that one can easily identify the Klingen 
Eisenstein series among a basis of eigenforms. And this is how we compute their Fourier 
expansions in practice. 

3.3. Maass Lifts 

This lift can be thought of as a modular form that can, in some sense, be computed easily 
from modular forms that are easier to compute. In particular, it can be derived from Jacobi 
forms and elliptic modular forms. The standard reference for Jacobi forms is [14) . 

Definition 3. A Jacobi form of level 1, weight k and index m is a function (p{z,T) for 
z eMi and r G C such that: 

(i) ^{^V7^)= + rf)'=e'^0(r, ^) for ( ^ ^) e SL(2, Z); 

(ii) (j){z, T + Xz + n) = e"2^""(^'^+2^^)0(z, r) for aU integers A, n; and 

(iii) (/) has a Fourier expansion 

n>0 r^<4nm 

A Maass lift |33j is a lifting from the space of Jacobi forms of weight k and index 1 to the 
space of Siegel modular forms of degree 2. 

Definition 4. Let </> be a Jacobi cusp form of weight k and index 1 and suppose (/)(r, z) ~ 
^jj^g (i(n, r)g"^'' is its Fourier expansion. Then the Maass lift of (j) is Maass(0) a Siegel 
modular form of degree 2 and weight k given by 

Maass(^)(-) = E ( E -^'"'^ ^) | ^''CV^ 

a,b,c \s\{a,b,c) ^ ^/ 

We observe that Maass((/)) corresponds to a formal Siegel modular form C and we compute 
C coefficient by coefficient: 

C^([a,6,c])= E '^'"'^ 

S \{a,b,c) 
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where d{n, r) are the coefRcients of the Jacobi form (f>. 

A map from ehiptic modular forms of level 1 to Jacobi forms can be found in jSOj ; we note 
that this is not the Saito-Kurokawa lift as it is not Hecke invariant. We summarize the result 
here. Let r] := J]^^^^(i - q") and 

A:=7?-6 J2 and 

r.s^Jj^r^s (mod 2) 
r,s£Z,r^s (mod 2) 



For / a degree 1, weight k modular form and g and degree 1 and weight fc + 2 cusp form, the 
map, / : (/, g) ^ |/A — {q-^f)B + gB is a Jacobi form of weight k and index 1. 

The composition Maass o / yields explicit formulas for the Fourier coefRcients of the forms in 
the Maass spezialschaar, i.e. in the image of the Maass lift, in terms of the Fourier coefficients 
of elliptic modular forms of level 1. The formulas for the Fourier coefficients of the Siegel 
Eisenstein series in Section IXTl are special cases of this construction. 



3.4. Theta Series 

Let g be a positive definite quadratic form on a lattice L of rank n and b{v, w) :— q{v + w) — 
q{v) — q{w) be its associated bilinear form. We assume q{L) C Z. The theta series of degree 2 
associated to L can be defined as 

e^^^Z):^ J2 e(tr(ZM[„,^])) 

{v,w)£L X L 

where M^y^^^ := Mf for the binary quadratic form / = q{xv + yw). It is well known that this 
series is a (scalar valued) Siegel modular form of degree 2 and weight n/2, for the level of the 
lattice and a quadratic character related to the determinant of the bilinear form. 

(2) 

The Fourier coefficients of 9)^ are related to representation numbers of binary quadratic 
forms by q. Namely, 

C{f) = # {{v, w)eLxL : / = q{xv + yw)} . (3.3) 

The standard properties of representation numbers imply that C{A ■ /) = C{f) for A E 
GL(2,Z). So, in particular, a theta series is a formal Siegel modular form with trivial formal 
character. 

The following algorithm computes the Fourier coefficients C(/) of a theta series of degree 2, 
for GL(2, Z)-reduced binary quadratic forms /. 



(2) 

Algorithm 1. Compute the Fourier coefficients C(/) of 9)^ for > disc/ > —X as well 
as C([0, 0,c]) for c < c„iax where c„iax is [^^J- 

1. [Precompute] For 1 < i < Cmaxi compute Vi := {v d L : q{v) = i} and Vi C Vi such 
that for every vector v € Vi, exactly one of w or — u is in V^. 

[Compute the singular coefficients] C([0, 0,0]) = 1 and for 1 < i < c^ax set C([0, 0, i]) — 

2-m , 

and for each a < 



2. 

3. [Iterate through all GL(2, Z)-reduced forms] For each 1 < a < 



c < 



4. 



4a 



look up Va and Vc- 

Iterate through vectors of shortest length] For v E Va, compute the linear form l{w) 

'2 if 6 7^0, 
4 if 6 = 0. 



b{v,w). For w g Vc, compute b = l{w) and X{b) 
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5. [Compute coefficient C[a, &, c])] If b<a and b^~4:ac>—X, then C[{a,b,c)] = 
C[(a,6,c)] + A(6). 

6. [Done] Return the map C{f). 

The algorithm is implemented in Sage a little differently (in particular the loops are in a 
different order in order to avoid computing things we do not really need to compute). 

3.5. Theta constants 

As we will see below in Section m being able to compute the products of theta constants 
is a very useful way to compute particular modular forms. For a, & € {0, 1}^ (column vectors), 
define the theta constant 9a,b as 

e^Az)^ E ^Q^M + J'^fe) (3.4) 

l=a mod 2 

where we use A[B] — *BAB. It is easy to see that 9a.b is identically zero if *a6 is odd. We note 
that an individual theta constant is not, in general, a formal Siegel modular form as can be 
seen from the following lemma: 

Lemma 3.1. For A e GL(2, Z), f eQ, we have: 

Ca,b{A ■ f) = CA-ia,tAb{f)- 

Note for example, that 0oooo(.^) and Oooio{Z) + 9oooi{Z) + Oqoh^Z) are formal Siegel 
modular forms but 0ooio(.^) is not. 

We have the following result that allows for the computation of the coefficients of a single 
theta constant. As noted above, this form is not necessarily a formal Siegel modular form. 

Proposition 3.2. Assume that *ab is even. Then the map Ca.b : Q — > Z given by 

i~a mod 2 

describes the coefficients of 9a,b{8Z). Here Q is as defined in Section \2.2\ a — *{ai,a2) is 
identified with the linear polynomial aiX + and the sum is over all linear forms £ in 
Z[X, Y] which are congruent to aiX + mod 2 and whose square equals f. 

The claimed formula becomes obvious if one writes £ = *(i?i,^2) and Z[£] = £\t + 2£i£2Z + 

Particular Siegel modular forms are often described as an algebraic combination of theta 
constants. The following proposition allows for the computation of the Fourier coefficients of 
such forms. 

Proposition 3.3. Let 

n d 

F{Z) :=E"^n^-...^».(^)' 

i=i j=i 
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where we assume that all ^aijbij are even. Then the map C : Q ^ Z associated to F{8Z) is 
given by 

n-l 

C(/) = ^a, (_i)(^i(''.i)+-+^<.(^-))/2. 

i—1 ii,...,i(i=aii....,aid mod 2 



The C in the previous two propositions can be formal Siegel modular forms. In that case we 
only compute the values of C at GL(2, Z)-reduced quadratic forms in Q, taking advantage of 
the extra structure a formal Siegel modular form has. We note, also, that whether or not C is 
a formal Siegel modular form has to be determined before the computations are carried out. 

The above definitions and propositions can be generalized to other levels. One should change 
these results in the expected way: e.g., the sums are mod N, the forms 9a^b,N{'^N'^Z) are formal 
Siegel modular forms, etc. In particular, we write 

9N,aAZ)^ E ^ (^^[^] + ' (3.5) 

l=a mod N 

where again we use A\B\ = *BAB. 



Remark 1 . A useful observation is that one can always take 6 = to get the space spanned 
by theta constants. In particular, the span of 9a,b.N{Z) is the same as the span of 9afl,N{Z)- 



3.6. Vector-valued theta series 

Following |25l we can define vector-valued theta series. Take a rank n even unimod- 
ular lattice i, let Pg C[X, Fjj be pluriharmonic (i.e., for every A£ GL(2,Z), P{AX) = 
det{A)''Symj{A)P{X)). Then we define 

9l p{Z) = E p (^ x-^ gTTi((x,x)T+{x,y)z+(y,y)T') ^ 

One can show that 9l,p £ Af„/2+fcj (r). After collecting terms, 9l,p is evidently a formal Siegel 
modular form because of the invariance of the pluriharmonic polynomials that make up the 
coefficients. 

To describe the generators of a certain space of vector-valued Siegel modular forms, it is 
useful to define the following two pluriharmonic polynomials. Let a, 6 € C" be such that their 
inner products (a, a) = (a, b) — (6, b) = (here the inner product of a and b is given by ^ a^foi, 
without complex conjugation). Consider the two polynomials 

k 

pW) ^ ^(2;,a)(^-'')(2/,a)''X'^F^-'' 

^il - t.ix,a)^^-^\y,ardet (^^^^^f^ 5)) ^ 

With this notation, 9^^ p3 G Md/2,j{T), O^ pU)^ ^ G Afd/2+fej(r). 

In order to compute these theta series, one needs to adjust Algorithm [T] appropriately. 



3.7. WronsA^ians 

We recall a useful construction in [2]. A Siegel modular form of degree 2 is a function in 
three variables, t, z and t'. Suppose we have 4 algebraically independent Siegel modular forms 
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Fi, F2, F3, F4 of weights fci, fc2 7 ^3, ^4, respectively. Then the function 



dFi 



dFi 
dr' 



k2F2 
dF2 



dF2 
dr' 



dFs 



kiFi 

dFi 
9r' 



(3.6) 



is a nonzero Siegel cusp form of weight /ci + A:2 + ^3 + ^4 + 3. See [2l Prop. 2.1] for more details 
and for the statement for arbitrary degree. 

Computationally speaking, computing the determinant {Fi, . . . ,^4} in the straightforward 
way is very expensive. Many products (36, computing the determinant in the naive way) of 
multivariate power series would need to be computed in order to find F . However, using basic 
properties of the determinant, one deduces 

Proposition 3.4. With notation as above, we have 

F = {F,,F2,F3,F4 = E {T^)Cf2 {T2)Cf, {T^)Cf, (T4) 

Ti = [ai,bi,ci]>0 
T2 = [a2,b2,C2]>0 

T3 = [a3,i)3,C3]>0 

T4 = ia4,64,C4]>0 





k2 


ks 


hi 


ai 


0.2 


03 


a4 


bi 


b2 


b3 


bi 


Cl 


C2 


C3 


C4 



eC^aiT + ^biZ + ^c,T'). (3.7) 



This is the basis for Algorithm [21 below. The subtle part in the algorithm is determining the 
bounds of the coefficients of the T;. These bounds have to be as tight as possible, since even a 
slight inaccuracy leads to subsequent iteration of many unneccessary Fourier indices. 



Proof of Proposition 



Equation (|3.7p follows from writing (|4.ip in the form 



F(Z)=det( E CFATi)(ll)eiaiT + biz + ciT') 



r2 = [a2,fc2.C2] 

and then using that the determinant function is multilinear. 



e{a2T + b2Z + C2t') 



□ 



The algorithm below to compute {-Fi, . . . ,^4} can be changed in the obvious way to work 
for arbitrary degree. 



Algorithm 2. Compute the Fourier coefficients C(/) of {Fi, . . . , F4} for > disc/ > -X 
as weU as C([0, 0, c]) where c < [^J . 

1. [Precompute quadratic forms and modular forms] Compute Qx := {/ G Q : > disc/ > 
-X}U{/= [0,0,c] G : c< L^J}; compute coefficients Cf,(/),CF,(/),CF3(/),CF4(/) 
for all f eQx- 
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2. [Compute the coefficient at /] For each 4-tuple (Ti, T2, Ta, T4) where Ti £ Qx and 
J:T^ = f. set 



C{f) = (^1)^^^ (^2)C^^3 (^3)^^^. (T4) X 

Ti = [ai,bi,ci]>0 

T2 = [a2,&2,C2]>0 

T3 = [a3,b3,C3]>0 
T4 = [a4,b4,C4]>0 

3. [Done] Return the map C{f). 



ki k2 fcs fc4 

fli 02 03 04 

&1 &2 &3 ^4 

Cl C2 C3 C4 



3.8. Rankin-Cohcn Brackets 

The attempt to construct a modular form as a Hnear combination of products of derivatives 
of two modular forms is known as the Rankin-Cohen construction. A complete investigation 
of this and related constructions has been given by Ibukiyama in [24] . See also [51 IllL I34j for 
characterizations of various Rankin-Cohen brackets. 

3.8.1. Satoh Bracket The Satoh bracket is a special case of this general construction. Satoh 
examined the case p = sym2 ® det''. Suppose F £ Mk{T) and G £ Mk'{T). Then 

[F,G]sym, = ^i^l^^z^-li^^zG) £ Mk+k'A^). (3-8) 

In order to implement the Satoh bracket in our package, we need first verify that the forms 
dzF and dzG are formal Siegel modular forms, even if they are not Siegel modular forms. 



Lemma 3.5. Let F be a modular form. Then dzF as defined above is a formal Siegel 
modular form. 



Proof Let Fz :— dzF and Cpz ■Q^'C[X,Y]j be such that CFz{[o-,b,c]) takes the 
coefficient of Fz at [a, 6, c] as its value. We need to show that Cpz is formal; i.e., that 
CFz{A-[a,b,c]) = {detA)>'A-CFz{[a,b,c]). 

A straightforward calculation shows that Cpz ([(a, b, c)]) = CpHa, 6, c]){aX^ + bXY + cY^). 
Now 

CpziA ■ [a,&,c]) = Cf(A- [a,b,c])A- {aX'^ + bXY + cY'^) 

= {detA)''CF{[a, b, c])A ■ {aX^ + bXY + cY"^) 
= {detA)''A ■ (CF([a, 6, c]){aX^ + bXY + cY^) 
= {det A)'' A -Cfz {[a, b,c]). 

□ 

Because of Lemma 13.51 we can compute the Satoh bracket as the difference of the product 
of two formal Siegel modular forms. 

Assuming that multiplication of vector valued Fourier expansions is implemented an 
algorithm can now be formulated as follows. 



Algorithm 3. 

1. For / and g compute the Fourier expansion with coefficients Cqzf{T) — C{T){X'^Tii 
XYTi2 + Y^T22). 
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2. Multiply the Fourier expansions according to formula p.8p . 



3.8.2. Ibukiyama's Brackets There are constructions found in [25 that are analogous to 
the Satoh bracket but whose images are in A/fc^4(r) and Alk.eiT). To give an indication of what 
they look like, we provide an example: let F E Mk{T) and G G AIk'{T). Then if k is even, 



k'{k' + 1) 



G X 



2^ 



\ 



' drdz 
d'^F I 2 d^F 
dz'^ ^ drdr 

2 



(fc' + l)(fc+l) 



dr'dz 
§Z£. I 



dF dG 
dF d8jJ^^dG_ dF 
OFdS^J^hGdl 

dr' dr' 



^'dFdG 
dz 



k{k + l) 



dz^ 



d^G 

2 d^G 
drdz 

2 ' 



drdr' 
, d'G 
"dr'dz 
d^G 
dP^ 



3.9. Borcherds Products 

In [5] , Borcherds introduced a multiplicative lift from whose image can be a Siegel modular 
form of degree 2. The Fourier expansions of these Borcherds products can be computed in 
polynomial time |18j . We refer the interested reader to the latter article, co-authored by the 
first author, for more details. 



3.10. Further remarks 

3.10.1. Restriction method We are going to discuss the approach that was used to obtain 
the results in [43j. In this work Siegel modular forms of degree 4 were considered, but in 
principle the algorithm can also be used to obtain Fourier expansions of forms for high level 
and degree 2. The idea is as follows. Suppose we restrict a Siegel modular form / : H2 — > C 
to the subspace : Hi —> EII2, r i-> sr + C for some positive definite symmetric s S M(2, Z) 
and some symmetric ( G M(2,Q). This will lead to an elliptic modular form whose Fourier 
expansion can be deduced from the Fourier expansion of /. Since many cusps for / o (/)^^ 
can be identified in IHI2, there are non trivial obstructions to the Fourier expansion of /. In 
[43j . where all cusps identified with each other, this provided enough restrictions to actually 
compute the initial Fourier expansions of a basis for a space of fixed weight. They could not 
give a proof that this method will work for all Fourier indices and all weights. In the case of 
level subgroups a further difficulty enters. There will be cusps that cannot be identified with 
each other in H2. 



3.10.2. Poor-Yuen Conjecture The Fourier-Jacobi expansion of a Siegel modular form can 
be used to describe a Fourier expansion. A formal Fourier-Jacobi expansion is a formal series 
J2 fmq"^ of Jacobi forms fm that have level m. Chris Poor and David Yuen conjectured that 
a formal Fourier-Jacobi expansion is the Fourier-Jacobi expansion a Siegel modular form for 
the full modular group if and only if the Fourier coefficients of the Jacobi forms satisfy the 
condition fm(n,r) — fn{m,r). It was recently announced by Poor and Yuen that have proved 
this result. 



3.10.3. Invariant method and Igusa's theorem on theta constants We fix a level N. For 
a, 5 e {0, . . . , — 1}^ (column vectors), define the theta constant 9M,a,b as in p.Sp . It is easy 
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to see that 6N,a,b is identically zero if a / and b / N have half integral entries and *ab is odd. 
[30j tells us that also the converse holds. There we can also find a proof that they are modular 
forms with respect to the following subgroup r{N,N) of Sp(2,Z): 

{7 = ( ^ g ) e Sp(2, Z) : 7 = I4 (mod N); A*B, C^D have diagonals divisible by 2N} . 

As we will see below in Section being able to compute the products of theta constants is 
a very useful way to describe modular forms. 
An important result is 

Theorem 3.6. Suppose that T{N,N) CT is a congruence subgroup. The graded ring of 
modular forms for T is the integral closure of the ring of those modular forms with respect to 
r that are formed as linear combinations of products of the theta series 0N,a,b- 

This method to compute rings of Siegel modular forms is impractical for large level. 
Nevertheless its connection to coding theory makes it interesting for general level. This can be 
found in [37] . 

3.10.4. The basis problem and theta series The basis problem, as formulated by Eichler 
[12 )|13j . asks whether the space of elliptic modular forms for fixed weight, level and character 
can be spanned by ^-series. The affirmative answer in the elliptic case partially generalizes 
to Siegel modular forms. The most complete treatment can be found in [3], where Bocherer 
proves the affirmative answer in the case of squarefree level and also gives a brief survey on 
the current effort in this area. In particular, he proves the following: 

Theorem 3.7. Assume that TV e N is squarefree and fc > 5. Let X be a character for To (N). 
Then all cusp forms in Mk{To{N), x) a,re linear combinations of 9-series. 

We also mention the very early treatment in |39| . There Ozeki proves that in the case of 
the full modular group one might restrict to unimodular lattices if and only if the weight is 
divisible by 4. 

There is a special case of this question. The critical weight for degree n forms is n / 2; in the 
case of degree 2 the critical weight is 1. It is expected that these weights behave particularly 
well with respect to the basis problem. In |54j it was proved that every Siegel modular form of 
degree 2 and weight 1 is a linear combination of O-sevies. At this point it might also be useful 
to the reader to mention the following results of the third named author. In (22) it was proved 
that there are no Siegel cusp forms of weight one for any group ro(iV) of Hecke type. 

Any modular form of weight less than the critical weight is called singular. In the case of 
degree 2 these are the half integral weight forms. Singular weights are well-understood for 
arbitrary degrees (see [17j V In particular, we know that all weight ^ forms associated to the 
full modular group only have nonvanishing Fourier coefficients for indefinite indices. Hence, by 
a theorem of Serre and Stark [49] they are given by ^-series. 



4. Rings of Siegel Modular Forms 

The problem of determining rings of modular forms has been of great interest. Complete 
descriptions of some rings are presented in Section 14.21 These results are almost entirely based 
on a combination of the theory of Borcherds products fsee 13. 9p . 0-series, and the invariant 
method that we discussed in Section l3.10.3l 
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Besides these results, there are methods that have been developed for cases where a set of 
generators has not yet been obtained. Although we focus on the descriptions mentioned above, 
we discuss some of these other methods in the next few subsections. Indeed, even for moderate 
levels, rings of Siegel modular forms become complicated. For applications it is sometimes 
sufficient to obtain a set of modular forms that spans a space of fixed weight. In these cases 
it is necessary to know the dimensions of these spaces. For this reason we also review known 
dimension formulas. 

4.1. Hecke eigenforms 

Hecke eigenforms are the modular forms that are of the greatest arithmetic interest. In order 
to compute eigenforms, one needs to compute a space of cusp forms, find the matrix for a Hecke 
operator and from that find the eigenforms. In this section we address issues surrounding this 
process. First, one needs to know when one has computed a basis for the space and to do this 
one needs to know the dimension of the space. Second, one needs to carry out the linear algebra 
computations described above. 

4.1.1. Dimension formulas In most cases only formulas for dimensions of spaces of cusp 
forms are available. This might already be sufficient, since many arithmetic questions can be 
reduced to this case. In Table [T] we give some dimensions of spaces of Siegel modular forms. 

4.1.2. Computing Hecke eigenforms Among all modular forms, the most arithmetically 
distinguished are the Hecke eigenforms. Fix a space of Siegel modular forms with basis 
{Fi, . . . , Fn}. Because the Hecke operators are a commuting family of linear operators, there 
is a basis {Gi, . . . , Gn} for the space made up entirely of simultaneous eigenforms. 

The forms Gi are determined computationally as follows. First, determine the matrix 
representation for the Hecke operator T2. Do this by computing the image under T2 of each 
basis element Fi . Build a matrix N that is invertible and whose jth row consists of coefficients 
of Fj at certain indices Qi,.. . ,Qn- To ensure that N is invertible we pick the indices one 
at time, making sure that each choice of index Qi increases the rank of the matrix. We then 
construct a matrix M whose jth row consists of coefficients of the image of Fj under T2 indexed 
by Qi, . . . , Qn- Then the matrix representation of T2 is M ■ N^^. 

Second, one finds the eigenspaces of the matrix for T2 and uses that to describe the basis 
Gi,...,G„. Typically only a single member of a given Galois orbit is described. See, for 
example, J51J. 




4.2. Particular Rings of scalar-valued Forms 

4.2.1. The full Siegel modular group In f27\ and [28], the five generators of the ring of 
Siegel modular forms and level 1 are identified: they are E4, Eq, xio, X12, X35 and are of weights 
4, 6, 10, 12 and 35. The first two are Eisenstein series and the last three are cusp forms. In [50j 
an indication of how to compute the four generators of even weight is given. One composes 




Mfc(Sp(2,Z)) [271150) 
Sk(^i(i)),N >3,k>5 |9l [38l[55) 



Ski^oip)), P prime , fc > 5 [201 

'S'4(ro(p)), p < 13 prime [41] 

S'3(ro(p)),P< 23 prime [41] 

S'2(ro(p)), p < 41 prime [41] 

SilroiN)), N >1 [22] 

Sfe j(Sp(2,Z)), j > 0, > 4 [53] 



Table 1 . Spaces and references for their dimension formulas 
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an explicit map from classical to Jacobi forms and the Maass lift from Jacobi forms to Siegel 
modular forms of degree 2. The implementation in [46] of these four generators follows [50j. 

In order to get the full ring of Siegel modular forms of level 1 without character, we require 
a method to compute X35- In |28j . it is given in terms of theta constants. In |19j . it is given as 
a Borcherds product. In [2], it is given as a Wronskian (see Section [377)) : 

X35 ^ -77r-^{E4,Ee,Xio,Xi2}- (4.1) 

Here E4, Eq, xio- X12 are normahzed so that Ce^HO, 0, 0]) = Ceb ([0, 0, 0]) = 1 and C^^j, ([1, 1, 1]) = 
C^,,([l,l,l]) = l. 

If one is content with just a basis for the vector space of Siegel modular forms of a fixed 
even weight, an observation by the first author is useful and would make, if true in general, the 
computation of such a basis significantly less expensive. In particular, in |45| he has observed 
that up to weight 172 one can span the space of weight k Siegel modular forms of degree 2 via 
forms that are the sum of products of no more than two Saito-Kurokawa lifts; see Section [3.31 
for more details on Saito-Kurokawa lifts. Since computing a Saito-Kurokawa lift is rather 
straightforward, multiplying a pair of Siegel modular forms is not terribly expensive, and 
adding many Siegel modular forms is very fast, this would represent a significant speed-up over 
what is done now. 

4.2.2. Congruence subgroups of small level In papers by Ibukiyama and his coauthors, the 
generators of some rings of modular forms of small level have been identified. Let Fq {£) be the 
Klingen subgroup of level £. For a subgroup F we use M, (r,x) for the ring of Siegel modular 
forms of degree 2 with respect to the group F and with character x (we omit the x if it is 
trivial). 

The groups whose rings of modular forms are known to be described in terms of generators 
are listed in Table [2] 



Group and Character 










Reference 


To (2) 










m 


r^(3) := {7ero(3) 




= 1} where V'3 ( c B ) 


= (di!j).f- -11 (^B) 


e To (3) 


m 


r^(4) := {7ero(4) 


1/14(7) 


= 1} where i/'4 (c B) 


= (d^).f- -11 (^B) 


e To (4) 


m 



Table 2. Subgroups F for which the generators of the ring M*(F) are known. 



4.2.3. Siegel modular forms of degree 2 and level 2 Consider the ring Af,(Fo(2)). In [2^ it 
is shown that the ring has 5 generators: 



X 
Y 
Z 

K 

Xl9 



94 

^2 



i,i))/4 



- ^^(0,0). (0,0) "^(0,0), (0,1) "^(0,0). (1,0) 
= (^(0,0),(0, 0)^(0, 0),(0,1)^'(0,0), (1,0)^(0,0), (1,1); 

= (^(0,1), (0,0) ~ ^('o,i),(i,o))^/16384 

= (^(o,i),(o,o)6'(o,i),(i,o)^'(i,o),(o,o)^'(i,o),(o,i)6'(i,i), (0,0)^(1,1), (1,1)) /4096 
{X,Y,Z,K} 



The generators are, respectively, of weights 2, 4, 4, 6 and 19. The first four are given in terms 
of theta constants and the last one is given as a Wronskian. 
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4.2.4. Siegel modular forms of degree 2 and level 3 The ring M, (Fq ^'^^ (3)) has 5 generators. 
Consider the quadratic forms 




and 



The modular forms of note are 



ai 
S3 

14{Z) 



9 A,, 

9e6 — 99 El , 

^ {c^ - d^)e(Xr{Q{x,y)Z)) where 



[xiy-i - xsj/i) + (x2y4 - Xiy2) and 
(a;i2/4 - x^yi) + (x3y2 - a;2y3) 



Xi4 := {ai,/33,(53,74}. 

We observe that these forms, respectively, have weights 1, 3, 3, 4, 14. With this notation, 
then, 

M,(r^^'(3)) =C[«i,/33,53,74]©Xi4C[ai,/33,53,74]. 

4.2.5. Siegel modular forms of degree 2 and level 4 Now, we consider r^''(4). First, as m 
[28] . let 



X5 ■— ^(O,O),(O,O)(-^)^(O,O),(O4)(-Z')0(O,O),(1,O)(-Z')0(O,O),(1,1)(2')0(O4),(O,O)(^) 



X e 



(0,1), (1,0) 



1,0), (0,0) 



1,0), (0,1) 



1,1), (0,0) 



1,1), (1,1) 



(Z). 



Then, as in [2], let 



fl/2 
fl 
92 

h2 

h 

Xii 



^(o,o),(o,o)(2.Z') 



12 



f(0,0),(0,0) + C'(0,1),(0,0) + "^(1,0), (0,0) 

^(0,0), (0,0) (2.2')'^ + 0(0,0), (0,1) (2^)"^ + 0(o,o),(i,o)(2.Z')'^ 
(6'(o,o),(o,i)(22')6'(o,o),(i,o)(2^)6'(o,o),(i,i)(2-Z'))^ 

{/l,.92,/l2,/3}- 



^(i,i),(o,o)(2^)'' 
?(o,o),(i,i)(2^)'' 



ThenM,(r^^(4))=C[/i,g2,/i2,/3,Xii]. 

4.3. Particular spaces of vector valued Forms 

The ring of all vector valued Siegel modular forms 0^. ^ Mkj{T) is not finitely generated. 
For this reason the vector valued weight j is usually fixed. The resulting module is finitely 
generated over 0^, Mfc(r). 
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Before we summarize the results on such modules of vector valued modular forms, we remark 
that vector valued Siegel modular forms are always cuspidal. This follows from a much more 
general result by Weissauer |54) . that covers Siegel modular forms of all degrees. 

We summarize the spaces that we can compute in Table [31 

4.3.1. Satoh's Theorem By Satoh |48) we know that 0feMfc_2(r) is generated by 6 
elements, all of which can be expressed in terms of Satoh brackets (see Section [3.8. ip . More 
precisely, he shows that 

Mfc,2(r) ^E4,Ee] ■ Affe-io(r) ® [E^^Xio] ■ Mk-uiT)® 

[E4,Xi2] ■ Mfc_i6(r)] ® [Ee,Xio] ■ C[Ee,Xio,Xi2]k-i6® 

[E6,Xl2\ ■ C[E(i,XlO, Xl2]fc-18 ® [XlO, X12] ■ C[X10, Xl2]fc-22- 

By C[Ai, . . . ,An]k we mean the module of weight k modular forms that can be expressed in 
terms of generators Ai, . . . , An- 

4.3.2. Ibukiyama's Theorems In |25l . generators for the rings Mk^A and Mkfi are given. The 
ring of forms of weight (fc, 4) is generated by 10 forms that are defined in terms of differential 
operators similar to the Satoh bracket mentioned above. The ring of forms of weight (fc, 6) are 
generated by a Klingen Eisenstein series, two theta series with pluriharnionics (see Section [3.6l 
for more information), and four forms that are defined via differential operators. 

5. A particular implementation 

In this section we describe our implementation of a package that can handle a wide variety 
of Siegel modular forms. In particular, one can multiply Satoh brackets by scalar- valued Siegel 
modular forms and one can multiply products of theta constants by other similar products. In 
short, our package allows for the multiplication of any two formal Siegel modular forms (even 
if the formal Siegel modular forms do not actually correspond to actual Siegel modular forms) . 
Of particular note is that our implementation handles Siegel modular forms compactly and 
efficiently. 

5.1. An implementation in Sage 

A Siegel modular form is implemented in the Sage package |46j as a formal Siegel modular 
form, i.e. as a map C : Q ^ R, where i? is a module (or ring) with a GL(2, Z)-action, such that 

A-\C{A.f) - C{f) 

for all A in GL(2,Z) and all / in Q. Such a map can be readily implemented via a 
class SiegelModularForms_class encapsulating a (Python-)dictionary whose keys are the 
GL(2, Z)-reduced quadratic forms below a certain bound. The methods of this class include 
multiplication, addition, Rankin-Cohen brackets, Hecke operators, etc. The bottlenecks for an 
effective implementation of these methods are multiplication and reduction of integral binary 
quadratic forms. 



Group and Weights 


Reference 


r, j = 2, fc 6 z 


m 


r, J = 4,6, A: e Z [25] 



Table 3. Spaces of vector valued Siegel modular forms whose generators are known. 
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Remark 2. We observe that our implementation of Siegel modular forms in Sage allows 
for the precision of a Siegel modular form F to be described in a number of ways. When we 
say F has discriminant precision X (here X is a positive integer), we mean that the keys of the 
dictionary mentioned above are the GL(2, Z)-reduced quadratic forms of discriminant greater 
than —X. 

There are a handful of other ways to describe the precision of a Siegel modular form (e.g., 
trace and box precisions) . In [46' we include a Siegel modular form precision class that handles 
translations between the various kinds of precision. 

For many theoretical purposes, it is best to describe the precision of a Siegel modular form 
in terms of dyadic trace [401 1411 l42l 1431 144j . Poor and Yuen have exploited this idea to prove, 
among other results, a nonvanishing theorem that can be thought of as a generalization of 
Sturm's bound for classical modular forms. 

5.2. Multiplication 

The multiplication of Siegel modular forms is expensive because it involves the evaluation 
of triple sums, whereas the reduction has to be applied in almost every call to a method of the 
class SiegelMod.ularForms_class since its instances only store coefficients for reduced forms. 

The obvious algorithm for the reduction is a variant of the Euclidean algorithm: applying 
repeatedly to a form / the two steps [a, b, c] i— [a, b%a, c] (where b%a denote the Euclidean 
remainder of b modulo a) and [a, b, c] i— >■ [c, 6, a] yields eventually the reduced form which is 
GL(2, Z)-equivalent to /. 

For a fast multiplication it is useful to provide separate implementations according to the 
types of the values of the formal Siegel modular forms. Operations on integer valued Siegel 
modular forms can clearly be handled faster than on those taking values, for example, in a 
polynomial ring over a number field. 

We point out that by "fast multiplication" we do not mean Karatsuba multiplication or the 
like. There are Karatsuba-like algorithms for multivariate power series rings, but they do not 
take advantage of the invariance of the coefficients for the Siegel modular form. The (potential) 
gains won by implementing a Siegel modular form as a multivariate power series rings that has 
a Karatsuba-like multiplication are outweighed by the great number more coefficients such an 
object would have to store and the number of redundant multiplications that would have to be 
done. There is currently no known Karatsuba-type algorithm that respects the invariance of 
the coefficients of a Siegel modular form, but this will investigated in future work by the first 
author. 

As a final remark we note that the set of all formal Siegel modular forms for a given coefficient 
ring i?, i.e. the set of formal Siegel modular forms C : Q — ^ ^(x) (x ^ linear character of 
GL(2,Z)), is naturally equipped with a multiplication. Namely, for two such formal Siegel 
modular forms Ci : Q ^ R{xi) ^nd C2 ■ Q ^ R{X2) the map C defined by 

C{f)= J2 CAh)C,{f2) 
/i,/2eQ 

/ = /l+/2 

defines a formal Siegel modular form C : Q R{xiX2)- Note that we would like to allow 
different characters since we would like to be able to multiply, for example, a formal Siegel 
modular form corresponding to even weight scalar-valued Siegel modular form (and hence 
with trivial character) with formal Siegel modular forms corresponding to scalar-valued Siegel 
modular forms of odd weight (and hence with character det). 

However, if xi is different from ^2 there is no natural sum of Ci and C2. On the other hand 
in our Sage implementation it is desirable to view formal Siegel modular forms as elements of 
an ambient algebra. This allows for an implementation which is consistent with Sage's internal 
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coercion system, which can take then over the necessary coercion steps to multiply (and add) 
formal Siegel modular forms with different coefficient rings if the coefhcients can be coerced to 
a common ring. 

To view formal Siegel modular of a given coefficient ring as elements of an algebra we proceed 
as follows. Let S be the group of linear characters of GL(2, Z). We view a formal Siegel modular 
form with values in i?(x) as a formal Siegel modular form with values in the group ring -R[S] 
equipped with the GL(2, Z)-action 

{A, ^ c(x)ex) ^ ^(^) ^•^(^) 

via the natural embedding R{x) ^ ^[■^1 which maps r to re^. Here runs through the natural 
basis of the group ring i?[S]. Formal Siegel modular forms with values in then naturally 
form an algebra over the base ring of R (which is either R if the action of GL(2,Z) on R is 
trivial or the base ring R' of i? if i? is a polynomial ring in two variables over R' equipped 
with the natural GL(2, Z)-action) via the Cauchy product mentioned above and the obvious 
addition. 

5.3. Data 

At |51| . one can find data for scalar and vector valued Siegel modular forms of level 1. In 
the future, similar data for Siegel modular forms of level greater than 1 will be posted. In 
particular, it is our goal to have data for all the rings and spaces listed in Tables [3] and [2] 

6. Generalizations 

A class similar to the SiegelModularForms_class would work also for other types of 
automorphic forms like Hilbert modular forms, Siegel modular forms of higher degree, 
orthogonal modular forms, etc. We formulate here in an abstract way the kind of object should 
be implemented to treat all kinds of higher rank automorphic forms in a unified way. 

Let G be a group, let i? be a module commutative ring and M be a monoid, both equipped 
with a G-action (g, r) i— > g.r respectively (g, m) i— >• g.m (such that, for each g in G the maps 
r I— > g.r and m M- g.m are homomorphisms of modules or rings and monoids, respectively). We 
set 

i?|Af|G --{C : M R : ^g eG,me M : C{g.m) = g.C{m)} 

Note that i?|M]'^ , for a ring i?, is a subring of the ring of power series i?|M]. 

In fact, the types of automorphic forms mentioned above and many more possess a Fourier 
expansion whose Fourier coefficients lie in _R|Af]'^ for suitable choices of M, R and G. Table|4] 
summarizes various examples. An implementation of these objects would again be a class built 
around a dictionary whose keys are representatives for the orbits in G\M . The effectiveness of 
such an implementation depends, of course, on a good reduction theory for finding distinguished 
representatives for the classes in G\M . Such a class would be highly desirable and encourage 
the implementation of more types of automorphic forms. An implementation for Sage by the 
first author will appear soon. 

Most useful for explicit computations are those graded (modules over) rings of automorphic 
forms of the types mentioned in the table where generators are known and can be computed 
effectively. There are known examples, in addition to the one of Siegel modular forms of degree 
2 already considered, of spaces of Hilbert modular forms, Jacobi forms, and Siegel modular 
forms of higher degree whose generators are known. More recent examples include Hermitian 
modular forms on the full hermitian modular group 

Tl := {M e A/(4, Zl) : M*JM = J} 
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Elliptic modular forms 



Group G 
Monoid M 
G-action on M 
Module (or ring) R 
G-action on R 



Vector-valued Siegel 
modular forms of 
degree n and weight 
k,j 



Group G 
Monoid M 

G-action on M 
Module (or ring) R 
G-action on R 



GL{n,Z) 

Set of semi-positive definite integral 
quadratic forms / in n variables 
(gj)^ f{{Xu...,Xn)g) 
¥[Xi,...,X„]j 

(g,p) ^ dctig)'' X pi{(Xi, ... ,X„)g) 



Hilbert modular forms 
of (parallel) weight k 
over a totally real 
number field L 



Group G 
Monoid M 

G-action on M 
Module (or ring) R 
G-action on R 



Set of totally positive or zero elements 
in the inverse different of L 
{9,0.) g'^a 

F 

(g,r) N{g)''r 



Hermitian modular 
forms over the 
imaginary quadratic 
field L 



Group G 

Monoid M 



G-action on M 
Module (or ring) R 
G-action on R 



GL(n,Zi,) 

Set of semi-positive definite integal 
hermitian forms / over L with n 
variables 

{g,f)^ f{(Xi,...,X,,)g) 

F 

(g,p) ^ det(g)'=p 



Jacobi forms of weight 


Group G 




k and index m S> (in 


Monoid M 


Set of {D,t) in -d'^ X i?"! s.t. D = 


where i9 denotes 




mod 4»™?~l, D = or -D > 


the different of L) over 


G-action on M 


{(e,x),{D,r)) ((£2_D,e(r-|-x))) 


a number field L 


Module (or ring) R 


F 


(cf. 0) 


G-action on R 


((e, x), r) H-> N{e)''r 



Table 4. Types of Fourier coeiEcients for various types of automorphic forms (F denotes a 

commutative ring) 



with J := (_°jq), and where Zl is either the maximal order in L = Q(-\/— T) [T], |10j or 
the maximal order in L = 3) IQ:. The generators can, in these cases, be obtained as 

Hermitian Maass lifts or Hermitian theta series. 
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